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The electronic properties ol boron-nitride nanoribbons (BNNRs) doped with a line of carbon 
atoms are investigated by using density functional calculations. Three different configurations are 
possible: the carbon atoms may replace a line of boron or nitrogen atoms or a line of alternating 
B and N atoms which results in very different electronic properties. We found that: i) the NCB 
arrangement is strongly polarized with a large dipole moment having an unexpected direction, ii) 
the BCB and NCN arrangement are non-polar with zero dipole moment, iii) the doping by a carbon 
line reduces the band gap independent of the local arrangement of boron and nitrogen around the 
carbon line, iv) an electric field parallel to the carbon line polarizes the BN sheet and is found to 
be sensitive to the presence of carbon dopants, and v) the energy gap between the highest occupied 
molecular orbital and the lowest unoccupied molecular orbital decreases linearly with increasing 
applied electric field directed parallel to the carbon line. We show that the polarization and energy 
gap of carbon doped BNNRs can be tuned by an electric field applied parallel along the carbon line. 



Single layer hexagonal boron-nitride (h-BN) 
nanosheets have been recently produced and placed 
on top of a Si02 substrate [l|. Several experimental 
groups have also fabricated free-standing h-BN single 
layers using sputtering of controlled energetic electron 
beams 0, y|. Unlike graphene, a h-BN sheet is a 
wide gap insulator, as is hexagonal BN, which is a 
promising material for opto-electronic technologies 
P, @, tunnel devices and field-effect transistors H, 
Small flakes of h-BN, i.e. BN nanoribbons (BNNRs), 
are semiconductors when hydrogen-passivated on one of 
the edges The band gap of BNNR can be tuned 

by using different types of atoms for edge passivation 
|ll| - [l4| . The zigzag BNNR has half- metallic properties 
when the boron edge is hydrogen-passivated, and the 
nitride edge is bare llj . Under a hydrogen-rich environ- 
ment, the zigzag BNNR with two-hydrogen-terminated 
edges becomes a ferromagnetic metal Moreover, the 
electronic and magnetic properties of the BNNRs could 
be modulated by applying an extra transversal electric 
field. Zhang et al used the local density approximation 
(LDA) to study the energy gap and showed that it can 
be tuned by applying a transverse electric field on a 
semi-infinite hydrogen passivated BNNRs with zigzag or 
armchair edges 0. Similarly by using the LDA, Park 
et al found that a transverse electric field decreases the 
band gap of armchair BNNRs monotonically while the 
band gap of zigzag BNNRs either increases or decreases 
depending on the direction and the strength of the 
applied field [T^l- Additionally, random doping BN with 
C atoms leads to spontaneous magnetization [isl . [l6j . 
Sai et al found piezoelectricity in a heteropolar nanotube 
depending on the chirality and diameter of the nanotube 
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which originates from the piezoelectric response of an 
isolated planar BN sheet |17| . 

In this study, we report first-principles calculations for 
the BNNRs which are doped by a zigzag line of carbon 
atoms in the middle and passivated by hydrogen atoms 
at the edges. We study different possible positions of 
the line of carbon atoms with respect to the BN-lattice 
and found that piezoelectricity and consequently the per- 
manent polarization in the BNNRs can be strongly en- 
hanced, tuned or even eliminated depending on the type 
of atoms (B or N) that surrounds the dopants. We found 
that doping BNNRs by carbon atoms decreases the band 
gap while the electric polarization of the doped BNNRs 
depends on the type of atoms (B or N) that surrounds the 
dopants. The NCB system is polarized in an unexpected 
way while BCB and NCN have zero dipole moments. We 
find that an external planar electric field applied along 
the C-line either reduces or increases the band gap and 
the dipole moment depending on its direction. 

Model and Method 

Our system is a finite-size h-BN sheet which is pas- 
sivated by hydrogens at its edges. A rectangular h-BN 
ribbon 3.1 x 3.4 nrn^ in size is taken as a perfect nanorib- 
bon, see Fig. [T] Total number of B and N atoms in the 
sheet is 418, in addition there are 56 H atoms at the 
edges. A zigzag line of 28 atoms is inserted in the mid- 
dle of the system replacing a line of h-BN atoms. We 
have studied three different positions of the carbon line 
(C-line) with respect to the BN-lattice as illustrated in 
Figs. [TJb-d) and compare the results with the undoped 
sheet (Fig. [Ija)). Three different ways of inserting the 
central C-line are: i) a C-line in the center with boron 
and nitrogen atoms on both sides (Fig. IHb)) which we 
name NCB, ii) a C-line in the center with boron atoms on 
both sides (Fig.[ljc)) which we name BCB, iii) a C-line in 
the center with nitrogen atoms on both sides (Fig. [ijd)) 
which we name NCN. Notice that all these systems are 
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passivated by hydrogen atoms at the edges to saturate 
the edge chemical bonds. In order to have equal sizes on 
both sides of the C-line, one has to simulate a non-exact 
rectangular flake, i.e. the right-up and down-left corners 
are not orthogonal, as illustrated by the two dashed rect- 
angles in Fig. [Ijb) having the same size. Therefore in 
the RHS (LHS) rectangles in NCB, BCB and NCN there 
arc additional B(N), B(B), N(N) atoms, respectively (sec 
Figs.HTc-d)). 

We employ density functional theory as implemented 
in GAUSSIAN (G09) [3 which is an electronic struc- 
ture packages that uses a basis set of Gaussian type of 
orbitals. For the exchange and correlation (XC) func- 
tional, the hybrid functional B3LYP is adopted in G09. 
For some particular cases, and for comparison purposes, 
we also used the LDA as implemented in the BigDFT [T9| 
package which uses real-space wavelets. In both cases, 
the self consistency loop iterates until the change in the 
total energy is less than 10~^ eV, and the geometries 
arc considered relaxed once the force on each nucleus is 
less than 50 me V/A. Using the 6-31G* basis set in G09, 
we expect that our calculation is capable to provide a 
reliable description of the electronic properties of the dif- 
ferent systems. Natural bonding orbital (NBO) analy- 
sis is also performed at the B3LYP/6-31G* level of 
theory 

Lattice mismatch 

We recall that in bulk li-BN the in-plane B-N bond 
length dRN = 1.446A has been determined by X-ray 
. Early ab-initio calculations [12] yielded 
I.35A. We obtain dsN =1.45A which agrees much 
better with experiment. On the other hand the bond 
lengths d]\[c and dsc between the inserted C-linc atoms 
and the N and B surrounding atoms, as well as the intra- 
chain bond length dec between C-atoms, depend on the 
specific configuration of the atoms in the NCB, BCB and 
NCN systems (see Fig.[T]). Wc summarized the results of 
our calculations in Table I. We found for the three doped 
BNNRs that dec < dBN < dsc- Notice that the rela- 
tive sheet bond lengths d^c =1.4lA and I.43A in NCB 
and NCN, respectively, are an indication of the partially 
double bond character of the N-C bonds at the central C- 
line is remarkable. The difference in bond lengths around 
the C-line in NCB are clearly visible in Fig. [ijb). The 
values dsc =1.55A and d^c =1.4lA lead to a lattice 
mismatch and symmetry breaking between the left and 
the right part of our system along the inserted C-linc. 
In Fig. [21 the arrows show the corresponding atomic dis- 
placement pattern for NCB. On the other hand the sys- 
tems BCB and NCN (Figs. [1] (c,d)) are symmetric with 
respect to the C-line. The changes in the bond-length 
and the dependence of the configuration on the inter- 
face leads to changes in the electronic polarization and 
consequently to a piezoelectricity effect. The results of 
this electro-mechanical coupling will be discussed in the 
following part of this section. 

Electronic charge distribution 

The perfect infinite BN sheet is a neutral sheet with 
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a uniform charge distribution far from the edges. In- 
deed, in an infinite two dimensional h-BN sheet the point 
group symmetry is D3/1 and there is no permanent dipole 
moment. However in finite BNNRs structure the sym- 
metry is broken which allows piezoelectricity, i.e. the 
appearance of a dipole under mechanical deformation. 
Therefore, the h-BN sheet is the most s imp le crystalline 
structure which allows piezoelectricity [2J^ . On the 
other hand, in BNNRs the finite extension of the sheet 
and the passivation of the boundaries by H atoms lead 
to the appearance of a permanent dipole moment. For 
the case of BN (see Fig. [TJa)) using B3LYP we obtain 

a dipole moment Pq = (—10.58, —1.98) D, (in this study 
the electric dipole moment vectors point from the nega- 
tive charge to the positive charge, e.g. see red arrow in 
Fig- Ha)), i.e. the system becomes piezoelectric. In the 
following wc investigate the influence of carbon doping 
on the polarization in BNNRs. 

Inserting of C-line alters the charge distribution mostly 
in the neighborhood of the C-line (see Table I), while fur- 
ther away the charge distribution approaches the one of 
a perfect BNNR. The different values of electronegativity 
for B(w 2.0), C(« 2.6) and iV(« 3.0) largely account for 
the charge redistribution and the concomitant formation 
of local dipoles. 

For the NCB the local dipoles for B-C and C- 
N bonds around the carbon line are directed to the 
right (Fig. [^b)) which polarized NCB. Surprisingly, al- 
though local dipoles are in the a;-direction, we found 
that the total dipole moment is mainly in the y-dircction 
(parallel to the carbon line), i.e. using B3LYP we find 

^ = (-9.37, -42.89)D. We attribute this effect to the 
presence of carbon atoms in the middle of Fig.[IJb). The 
C atoms which are bonded to the B (N) atoms absorb 
(give) electronic charge yielding two different group of C 
atoms, i.e. those bonded to B atoms with negative elec- 
tronic charge and those bonded to N atoms with positive 
electronic charge. This is the main reason for the unusual 
direction of the dipole moment. In fact the magnitude of 
Py is related to the carbon line. Notice that this is seen 
only when different types of atoms surround the C-line, 
i.e. the BCB and NCN systems have small Px and Py. 

As seen from Figs.[IJc,d) in the BCB and NCN systems 
the dipoles of the left and right hand sides of the C-line 
are in opposite directions and no overall dipole moment 
in the x-direction arises from these edges. There is a 
small dipole moment Py because the corners are differ- 
ently terminated at the left and right edges. 

A similar description is possible for the strong local 
dipoles on the B-C and N-C bonds induced by the inser- 
tion of the C-linc. Although having a C-line in the mid- 
dle, the BCB and NCN systems have vanishing dipole 
moments in the x-direction because the local dipoles on 
the B-C and N-C bonds around the C-line are directed 
oppositely. The direction of the local dipole moments, 
indicated by the horizontal black arrows in Figs. [Ijb-d), 
is chosen according to the electronegativity differences 
of the atoms as discussed earlier. It is important to note 
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that the number and arrangement of hydrogen atoms are 
the same in all model systems, thus the main reason of 
the strong total polarization is the local arrangement of 
B and N atoms around the C-line and not the presence 
of the H-atoms of the edges. 
Electrostatic potential 

The electrostatic potential (ESP) mapped on the plane 
of the h-BN sheet for the representative models is shown 
in Fig. [3l The sign of ESP at any point depends on 
whether the ESP due to the nuclei and the electrons is 
dominant at that point. The difference in electronegativ- 
ity between the B and N atoms results in a combination 
of weak ionic and covalent bonding in the h-BN sheet. 
Fig. ^a) shows that the B-N bonds have a larger ESP 
(blue spots) around the B ions as compared to the N 
sites (red spots). It shows that the N ions attract elec- 
tronic charge from the B ions as one expects in view of 
the atomic electronegativities. A characteristic pattern 
of alternating positive and negative ESP regions is seen 
over the BN surface in which the left side has generally 
less ESP than the right side. Accordingly, this explains 
why the BN has a permanent dipole moment in accor- 
dance of the electronegativity differences on H-B and H- 
N bonds at the edge which saturate the dangling bonds. 
The NCB, as shown in Fig.[3l^b), has a steep gradient in 
the ESP near the C-line which results in a strong dipole in 
this region which confirms the direction of the dipole mo- 
ment shown in Fig. [ijb). The big arrows in Figs. |3Ia,b) 
refer to the dipole moments. Similarly, the direction of 
the opposite dipolcs in Figs.[TJc,d) can be well explained 
by the ESP variation in Figs. [3l(c,d). Strong variation of 
the ESP around the carbon lines causes strong dipoles 
but directed oppositely. In NCB, the antisymmetrical 
distribution is also due to the chosen non-perfect geom- 
etry of the corners. 

Polarizability 

The polarization of the BNNR system varies in re- 
sponse to an external electric field. Here we focus on the 
NCB system which has an appreciable permanent polar- 
ization. We apply a uniform electric field in the plane 
of the flakes parallel to the C-line (y-direction) and al- 
low them to relax each time the electric field is altered. 
The reason for directing the field along y is that we want 
to compensate the largest component of the spontaneous 
dipole moment Pq, i.e. Py. The field strength E = Ey 
is then increased step-wise until the induced dipole mo- 
ment completely compensates Py. We denote the com- 
pensating field strength by Eq. Wc also inverted Ey to 
check the linear dependence. Red symbols in Fig. U (for 
NCB system) show that the dipole moment has almost 
a linear dependence on E within the studied range (- 
1,0.8) V/nm. The slope of this line gives the polarizabil- 
ity ayy = dE/dPy. Yiovti B3LYP we found ayy—^b.Q& D 
nm/V and i^o ~ 1 V/nm. In Figs. |H[b,c) we show the 
ESP and Pq for £'^=-1.0 V/nm and £^=-^0.82 V/nm, re- 
spectively. The lattice deformations due to applied elec- 
tric field are shown in Figs. |ll[d,e). It is interesting to 
note that a positive electric field shifts the RHS-atoms 



to the bottom and the LHS-atoms to the top, and vise 
versa for negative field. This yields a shear stress along 
the C-line (filled black symbols). The larger the electric 
field the larger deformation and the larger the enhance- 
ment of piezoelectricity. 
Metallic Properties 

To show the effect of the insertion of the C-line on 
the conductivity of BNNRs, we have listed the calcu- 
lated energy gaps for our systems in Table H] The gap is 
calculated as the energy difference between the highest 
occupied molecular orbital (HOMO) and the lowest unoc- 
cupied molecular orbital (LUMO). It is well known that 
LDA (which we used in our BigDFT calculations) under- 
estimates the energy of the excited levels and therefore 
the HOMO-LUMO gap. The hybrid functional B3LYP 
is expected to give a better estimation. However, both 
approximations reveal the same trend i.e. significant de- 
crease in the HOMO-LUMO gap after insertion of the C- 
line. It implies that the electric conductivity is increased 
in the doped sheets, however their gaps are still large 
and therefore the system is not metallic. We will show 
later that the HOMO-LUMO gap can also be tuned by 
an electric field. Figure [5] illustrates the density of states 
(DOS) for the studied BNNRs according to B3LYP. Note 
that the insertion of a C-line generates a few new elec- 
tronic states within the wide gap between the HUMO 
and LUMO of the perfect BN, reducing the gap width 
significantly. Notice that the distribution of energy levels 
around the HOMO-LUMO levels in BOB and NCN are 
different and the energy levels for NCN (BCB) are sig- 
nificantly separated below (above) the HOMO (LUMO). 

The change in the metallic properties can be explained 
by looking at the distribution of the electronic states near 
the Fermi level. As seen from Fig. |3] in the BN system 
the HUMO (LUMO) is localized at two corners around 
the N (B) atoms. In the BN system, the HOMO is lo- 
calized at the right-down side with zig-zag edges while 
the LUMO is localized at the left-upper corner with zig- 
zag edges where these corners are orthogonal opposite to 
the other corners, i.e. right-up and down-left with zig-zag 
edges. It shows that both chirality and the local geometry 
of the edges arc important parameters for localizing the 
electrons. On the other hand in the doped systems the 
HOMO and LUMO states are localized around the C-line 
over either B-C or N-C bonds, see Figs. [3i;b,c). In BCB 
and NCN systems the C-line opens a wide channel of few 
atoms through which the electronic charge can be trans- 
mitted. But the conductivity in the direction perpen- 
dicular to the C-line is expected to be small because no 
contribution from HOMO and LUMO is extended along 
this direction. Notice that the localized orbitals in BCB 
and NCN have different patterns. In BCB both HOMO 
and LUMO are localized on the C atoms and the adja- 
cent B atoms (i.e. on the B-C bonds) , while the N atoms 
are not contributing. This is apparently in contrast to 
the belief that the HOMO (LUMO) are locafized on N 
(B). This is due to the larger electronegativity of C as 
compared to B. In the NCN system HOMO and LUMO 
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are localized mainly on the C-C bonds (with different 
orientation for HOMO and LUMO) and are partially lo- 
calized on the adjacent N atoms and the B atoms never 
contribute. 

The special shape of the localized HOMO and LUMO 
states around the two ends of the C-line for the NCB is 
of most interest. The HOMO orbitals in NCB arc lo- 
calized over a few top C-C bonds and neighbor B atoms 
while for the LUMO we notice that they are localized 
over a few bottom C-C bonds where B (N) atoms in the 
RHS (LHS partially) are mainly contributing with differ- 
ent shape as compared to the HOMO. It is important to 
emphasize that: i) we found the same orbital distribu- 
tion in a BN system when it is doped by an arm-chair 
line of carbon atoms in the middle (those results are not 
reported in this paper), and ii) we did not find such or- 
bitals employing semi-infinite NCB, i.e. a ribbon with 
periodic boundary along the C-line. The special shape of 
HOMO and LUMO orbitlas found in NCB are not found 
in the BN system 

Next we will focus on the NCB system to study its 
response to an external electric field. When an electric 
field in the opposite direction of the dipole moment of the 
NCB is applied, the gap becomes smaller, see Fig. El^a), 
and the HOMO (LUMO) is localized much closer to the 
top (bottom), see Fig.[nic). On the other hand, we found 
in the previous section that the electric polarizability in 
response to an external field is increased by doping the 
BNNRs by carbon atoms. This trend can be used to 
obtain metallic properties in doped NCB. Insertion of 
different numbers of carbon lines as well as using differ- 
ent arrangements of the atoms near the C-line can be a 
practical tool to tune the electric transmission through 
the BNNRs and their metallic properties. 

An external uniform electric field Ey along the C-line 
in NCB causes a potential difference along the C-line. 
We assume that the HOMO and the LUMO are effec- 
tively centered around yn and respectively An ap- 
plied electric field (Ey) can modulate the energy gap 
as (j/L — yH)sEy, where e is the elementary charge. 
The gap can be tuned by the strength and direction 
of the electric field. The latter linear behavior is ob- 
tained even using the LDA (see top inset in Fig. IDJa)) 
and is continued (with a slope of (j/i — yi/)e) while the 
centers of the HOMO and the LUMO are more or less 
fixed which is the case even for a rather strong field of 
E=+0.82V/nm (compare Fig. 3(b) and Fig. 6(c)). The 



same linear behavior is found for the reverse electric fields 
Ey > -0.8V/nm, see Fig. 6(a). 

On the other hand, if a strong negative field about 
— IV/nm is applied, the dependence is no longer lin- 
ear. Fig. 6(b) shows that such a negative field alters 
the spatial distribution of the orbitals and elongates both 
HOMO and LUMO into a non-localized distribution sim- 
ilar to those found for the BCB or NCN systems. There- 
fore around the maximum, ~ yn and the energy 
gap is independent of the field strength. However for 
the negative fields beyond the maximum the energy gap 
decreases which corresponds to the localization of the 
LUMO and the HOMO where yL < yn- This simple 
picture argues that {yL — yH)^ Ey modulates the energy- 
gap. Notice that the slope {yL — yH)e changes sign pass- 
ing the maximum {Ey ~ -1 V/nm), therefore by decreas- 
ing the electric field, the gap decreases when Ey > - 
IV/nm and increases when Ey < -1.0 V/nm. A similar 
behavior was found for a semi-infinite non-doped BN sys- 
tems 0, , with difference in maximum and energy gap 
values. The bottom inset in Fig. 6(a) shows the vari- 
ation of the HOMO and LUMO sates with the electric 
field. Wc believe that the control of the conductivity by 
an external field as wc found for the NCB will be of great 
practical interest. 

In summary, we studied the electronic properties of 
boron-nitride sheets in the presence of a zigzag atomic 
line of carbon atoms in the middle of the BN sheets. 
Electronic polarization of the BNNRs was found to 
depend on the local arrangement of the boron and 
nitrogen atoms around the carbon line. Doping with 
carbon atoms decreases the non-metallic properties of 
the BNNRs by reducing the energy gap and increasing 
the polarizability. By applying an electric field along 
the carbon line we showed that the dipole moment of 
BNNR induced by the carbon dopants can be reduced 
and even eliminated. This also reduces the energy gap 
in the electron spectrum and allows one to control the 
conductivity with an external field. 
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TABLE I: Bond lengths in A, atomic charge in e, HOMO-LUMO band gap in eV and spontaneous electric dipole moment Pq in 
D for the four studied nanoribbons in Fig. [1] Results for the gap are given for both B3LYP and LDA (respectively calculated 
by G09 and BigDFT). 
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FIG. 1: (Color online) (a) Perfect boron-nitride nanoribbon (BN). (b) Boron-nitride nanoribbon with a zigzag C-line in the 
middle surrounded by both boron and nitrogen atoms on both sides (NCB). (c) Boron-nitride nanoribbons with the C-line 
surrounded only by boron atoms on both sides (BCB). (d) Boron-nitride nanoribbons with the C-line surrounded only by 
nitrogen atoms on both sides (NCN). In (a) and (b) red arrows refer to the total dipole moment of the systems and in (c), (d) 
red arrows refer to the total dipole moment on the two sides. The black horizontal arrows refer to the local dipoles corresponding 
to the B-C or N-C bonds. The simulated systems has size 3.46x 3.11 nm? (Fig. [l]). Numbers refer to the bond lengths in 
A unit. 
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NCB 

FIG. 2: (Color online) Atomic displacement after inserting an atomic carbon line in the middle of the BN sheet, i.e. NCB. The 
color coding indicates the y-component of the displacement vector, i.e. red (blue)= downward (upward) displacement. 
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(a) BN: ESP 



LUMO 



HOMO 




(b) NCB: ESP 



LUMO 
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(d) NCN: ESP 
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FIG. 3: (Color online) Contour plots of the electrostatic potential (ESP) for the four systems of Fig. [T] and the corresponding 
highest occupied molecular orbitals (HOMO) and the lowest unoccupied molecular orbital (LUMO). In (a) and (b) the black 
arrows refer to the total dipole moment. 
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(a) 




. 
x(A) 

FIG. 4: (Color online) (a) Electric dipole moment versus applied electric field, both in the y-direction, for the NCB system 
using B3LYP functional. The corresponding ESPs and lattice deformation for E — Ey=-1.0 V/nm (b,d) and E = £^=+0.82 
V/nm (c,e). The black arrows (symbols) refer to the direction of the applied electric field (carbon atoms) and the white arrows 
refer to the total dipole moment. 
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FIG. 5: (Color online) Density of states (DOS) for the four model systems. The presence of the carbon lines in the middle of 
the systems decreases the gap and modifies the DOS in different ways for each particular cases. 
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FIG. 6: (Color online) (a) The HOMO-LUMO energy gap for NCB versus applied electric field using B3LYP. Bottom inset 
shows the change in HOMO and LUMO versus applied electric field and the top inset shows the LDA results for positive electric 
field. The change in the HOMO and LUMO orbitals of NCB in the presence of electric field Ey =-1.0 V/nm (b) (+0.82 V/nm 
(c)). 



